Accepted to ApJ 

SLIP-SQUASHING FACTORS AS A MEASURE OF 
THREE-DIMENSIONAL MAGNETIC RECONNECTION 

V. S. Titov^ 

Science Applications International Corporation, 10260 Campus Point Drive, San Diego, 

California 92121-1578 

titovvQpredsci . com 

and 
T. G. Forbes 

University of New Hampshire, EOS Institute 8 College Road, Durham, NH 03824 

terry . f orbesOunh . edu 
and 
E. R. Priest 

School of Mathematics and Statistics, University of St. Andrews, Fife KY16 OSS, Scotland, 

UK 

ericSmcs . st-and .ac.uk 

and 

Z. Mikic\ J. A. Linker^ 

Science Applications International Corporation, 10260 Campus Point Drive, San Diego, 

California 92121-1578 

mikiczOpredsci . com, linker jOpredsci . com 
ABSTRACT 

iJn transition to Predictive Science, Inc., 9990 Mesa Rim Road, Suite 170, San Diego, CA 92121-2910 



-2- 



A general method for describing magnetic reconnection in arbitrary three- 
dimensional magnetic configurations is proposed. The method is based on the 
field-line mapping technique previously used only for the analysis of magnetic 
structure at a given time. This technique is extended here so as to analyze the 
evolution of magnetic structure. Such a generalization is made with the help 
of new dimensionless quantities called "slip-squashing factors". Their large val- 
ues define the surfaces that border the reconnected or to-be-reconnected magnetic 
fiux tubes for a given period of time during the magnetic evolution. The proposed 
method is universal, since it assumes only that the time sequence of evolving mag- 
netic field and the tangential boundary fiows are known. The application of the 
method is illustrated for simple examples, one of which was considered previously 
by Hesse and coworkers in the framework of the general magnetic reconnection 
theory. The examples help us to compare these two approaches; they reveal also 
that, just as for magnetic null points, hyperbolic and cusp minimum points of 
a magnetic field may serve as favorable sites for magnetic reconnection. The 
new method admits a straightforward numerical implementation and provides a 
powerful tool for the diagnostics of magnetic reconnection in numerical models 
of solar-flare-like phenomena in space and laboratory plasmas. 

Subject headings: MHD — Sun: corona — Sun: magnetic fields 



1. INTRODUCTION 

In spite of the progress reached so far in understanding of the key role of magnetic re 
connection in many phenomena of space and laboratory plasmas, its quantit ative description 



in re a listic three-dirn ensional (3D) configurations still remains a challenge (jPriest fc Forbes 



2000l : ISchindle]l2006h . This concerns, in particular, the global aspects of reconnection such 



as an identification of reconnected magnetic fiuxes and related structural changes in solar 
coronal configurations. There are basically two fundamental reasons why this important 
problem is so complicated in the 3D case. 

First, by analogy with the two-dimensional (2D) case, one could assume that partition- 
ing the configurations into distinct interacting magnetic fiuxes by separatrix surfaces (SSs) 
and estimating their time variation might help to resolve this problem. And, indeed, this 
idea works very well whenever such SSs are due to the presence of magnetic null points (NPs) 
in the solar corona — the NP SSs are boundless surfaces, and so they provide the required 
partition. However, the situation is different if the SSs are due to the presence of so-called 
bald patches (BPs), which are the segments of the photospheric polarity inversion line where 
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coronal field lines touch the photosphere (|Seehaferl Il986l : iTitov et al.l Il993l : iBungey et al. 
19961 : iDelannee fc Aulanierlll999l : lAulanier et al.ll2002l ). The BP SSs are usually not bound- 
less and, therefore, not able to provide the required partition of the respective configura- 
tions. It is also important that they are not exceptional structures, but are highly likely 
to be present in practice in the comp l ex configurations that exist in the solar atmosphere 



(ITitov et al. 



; m p 
I993I: 



Bungev et al. 



1996 



Titov fc Demoulin [1999; Delannee fc Aulanier 1999 



Aulanier et al. 



2002 : 



Aulanier fc Schmiedei 



2002 



Pariat et al.ll2004l ). Moreover, the set of 



3D structural features where reconnection may occu r is even wider than mentioned above, be- 
cause it includes also quasi-separatrix layers (QSLs) ( Demoulin et al.|[l996 : Priest fc Demoulin 
19951 ). Similar to BP SSs, these QSLs also do not provide an exact partition of configurations 
into distinct magnetic fluxes. 

The second fundamental difficulty in quantitative description of 3D reconnection is 
related to magnetic diffusion and the breakdown of frozen-fiux behavior. Being an inherent 
part of the reconnection itself, magnetic diffusion in 3D can occur throughout the plasma 
volume in evolving configurations. This is of particular importance for numerical MHD 
simulations, which may presently be performed only with a relatively high numerical or 
physical diffusion compared to real solar conditions. Deviations from the frozen-in law 
condition are accumulated in this case throughout the volume to form a background against 
which it becomes difficult to distinguish between magnetic diffusion and reconnection. 

Nevertheless, significant progress in understanding 3D reconnect ion has been achieved 



in the framework of general magnetic reconnection (GMR) theory (jSchindler et al.l Il988 



Hesse fc Schindleiill988l : iHesse et al.lll99ll : iHesse &: Birnlll993l ). which is focused on the anal- 
ysis of magnetic field ev olution constrained only by Maxwell's equations and Ohm's law. On 
the basis of this theory, iHesse et al.l (120051 ) have recently shown that, regardless of the pres- 
ence of certain structural magnetic features, the time change of the reconnected magnetic 
flux is always directly related to the maximum of the field-line integral of parallel electric 
field. For a plasma-magnetic configuration contained in a closed volume, this means that the 
change of the magnetic connection between the boundary plasma elements is fully controlled 
by the voltage drop along these lines. It is clear, however, that such a change of the magnetic 
connection can always be found without calculating the electric field if the magnetic field in 
the volume and the plasma flows on the boundary are known at every moment. As will be 
shown, the theory based on this premise complements the GMR theory by providing new 
intriguing opportunities for the analysis of magnetic reconnection. The aim of this paper 
is to develop a general mathematical formalism that allows one to describe the change of 
magnetic connectivity and identify the reconnecting r nagnetic fluxes in any evolving config- 
uration. The new theory has been outlined earlier by iTitovl (l2007bl ). This paper, however, 
provides its full exposition. 
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The new formalism can also be viewed as a generalization of the theory of coronal 
magnetic connectivity, whose spatial variation is measured with the help of the so-called 
squashing factor Q. This quantity wa s introduced first for closed magnetic configurations 



with a plane photospheric boundary (ITitov et al.l Il999l : iTitov &: Hornigl l2002l : iTitov et al. 



20021) and th en extended to other configurations with arbitrary shapes of the boundaries 
( lTitovll2007al ). The Q factor allows one to identify in a given configuration at a fixed time 
all the structural magnetic features, such as SSs and QSLs as well as their hybrids. Here we 
generalize this quantity in order to characterize the evolution of magnetic connectivity by 
applying the Q factor to so-called slip mapping — the mapping that describes the slippage 
between boundary plasma elements and footpoints of magnetic field lines. 



It was previously noticed by lForbes &: PriestI (119841 ) for 2D magnetic configurations that 
the speed with which the photosphere is swept by a magnetic separatr ix line allows one to 
estimate the electric field at the reconnection site. iForbes fc Liru (120001 ) have also proposed 
strong arguments that this principle can be extended to 3D as well, excluding the above- 
mentioned configurations with QSLs, if one focuses only on the amount of magnetic flux 
reconnected up to a particular instant of time. Our paper provides exactly this type of 
generalization of the principle for 3D configurations with QSLs included. 

Section [2] describes the general theory of evolution of magnetic connectivi ty, which is 



i llustr ated in section [3] for the example of magnetic evolution proposed earlier by lHesse et al. 



( 120051 ). Section m presents another example of magnetic evolution in a current layer config- 
uration to show the difference between reconnection and magnetic diffusion. In section El 
we discuss the general properties of magnetic minimum points, which happen to be closely 
related to the considered theory. The obtained results are summarized in section [61 



2. EVOLUTION OF MAGNETIC CONNECTIVITY 

The basic concept that makes it possible to define the frozen-in law and its violation 
i n magnetized -plasma flows is the connection of plasma elements by magnetic field lines 



( lAxfordll 19841 ). When trying to characterize the magnetic structure and its evolution, we find 
it extremely useful to apply this concept to the plasma elements that are currently located at 
the boundary of a given configuration. The magnetic field lines threading the plasma volume 
generally start and end up at some pieces of the boundary with the coordinate charts {u^, v?) 
and (w^,w^) (Figure ill). In this way, they define a mapping (m^,m^) {w'^^w'^) determined 
by some vector function {W^ [u^ , u^) ^W"^ {u^ , u^)) from one piece of the boundary to the 
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other. The local properties of this mapping are described by the Jacobian matrix 



D 



(1) 



For each field line, this matrix determines a linear mapping from the tangent plane at the 
launch footpoint to the tangent plane at the target footpoint, so that a circle in the first 
plane is mapped into an ellipse in the second plane (Figured]). The aspect ratio A of such an 
ellipse defines the degree of local squashing of elemental fiux tubes, which means that any 
infinitesimal circle centered at a given launch point is mapped along the field lines into an 
infinitesimal ellipse with this aspect ratio at the target footpoint. 



Fig. 1. — Circle is mapped into an ellipse by a linearized field-line mapping acting between 
the tangent planes of the launch and target boundaries, where two different curvilinear 
coordinates and respectively, exist. The aspect ratio of the ellipse, when it 

is large, coincides with a high value of the squashing factor Q. 

The value A is invariant with respect to the direction of the field-line mapping and so 
it can be used for characterizing the spatial variation of magnetic connectivity. It turns out, 
however, that Q = X + A^^ provides a more compact expression in terms of coordinates 
and practically coincides with A in the most interesting cases where it is large. So it is 
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natural to use the value Q, called the squashing degree (ITitov et al.l 12002) or squashing 
factor, for characterizing the local properties of magnetic connecti on between the boundary 



elements. The Q fact or can be expressed in terms of the norm N (jPriest fc Demoulirull995 



Demoulin et al.lll996l ) and Jacobian A of the field-line mapping as follows 

Q = ivV|A|. 



(2) 



It is con venient for fur ther consideration to rewrite the covariant expressions for and A 
found by iTitovl (l2007al ) for an arbitrary shape of the boundary in the following matrix form 



iV2 

lAI 



tr {D^G*DG-^) , 
I det {D'^G^DG-^) 







Bn 






B* 

n 



(3) 
(4) 



where Bn and B* are the normal components of the magnetic field at the conjugate launch 
and target footpoints, respectively, and 



G 



dR OR 



(5) 



is the matrix of the covariant metric at the launch footpoints with det G denoted by g. The 
radius-vector R{u^,u'^) defining the shape of the boundary in a chosen coordinate chart 
(u^, u^) is assumed here to be known. The latter equally refers to the boundary at the target 
footpoints with radius- vector where the matrix of the covariant metric is 

'OR dR' 



G* 



(6) 



The asterisk symbol (superscript) indicates here that the values of G refer to the conjugate 
target footpoints such that {w^^ w'^) = (W^{u^, v?), W'^{u^ .u^)). In other words, the asterisk 
symbol pulls back the values of G from the target to the launch footpoints at which the Q 
factor is evaluated. 

If the magnetic field B(r) is defined analytically or on some numerical grid, the Jacobian 
matrix ([1]) can be straightforwardly computed from the coordinates of the target footpoints 
that are obtained by integrating the field-line equation 



dr 

d7 

with the initial condition at launch footpoints 



B 



(7) 



It=0 



Riu^. 



u 
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and the stopping condition at target footpoints 

rl^^^ = Riw\w'), (9) 

which are satisfied by some parameter values Tg to be found. Tracing of the field lines 
should be done with a sufficiently high accuracy to make the subsequent calculation of the 
coordinate derivatives accurate enough. Although such computations are relatively time 
consuming, they may be fulfille d with a very goo d resolution even for the entire solar corona 



without using supercomputers (ITitov et al.ll2008l ). The latter are required only for numerical 



modeling of the coronal magnetic field itself. For example, such a computation for the given 
magnetic field data on a Mac Pro with one processor Intel Xeon (2GHz) requires only 23 
hr. The resulting Q factor is obtained on a uniform grid with 3600 and 7200 mesh points 
in latitude and longitude, respectively, for a magnetic field defined in spherical coordinates 
(r, 6, (f)) on a nonuniform grid with 174 x 118 x 174 mesh points. 

As mentioned in section [T|, the Q factor allows one to identify in a given configuration 
all topological and geometrical magnetic features at a fixed time instant. Formally, the SSs 
are distinguished from the QSLs by infinite rather than large finite values of Q. However, 
from a numerical point of view, these features are similar, since they are both identified via 
large finite values of Q. The only difference is that the maxima of Q increase at SSs limit- 
lessly with a decreasing size of a numerical grid used for computing Q, while at QSLs such 
maxima converge to certain large values. Such a subtle difference is usually not important 
for applications and so, for brevity, we will use hereafter the term QSLs in a broad sense to 
include also SSs, unless it is stated otherwise. 



2.1. Ideal evolution of the magnetic field 

If the evolution of magnetic configurations is subject to the frozen-in law, the initial 
connection of plasma elements by magnetic field lines remains unchanged. The analysis of 
magnetic connectivity then is significantly simplified, because the field-line mapping Ilf at 
any time t can be composed in this case as follows: 

= FionooF„i, (10) 

where the superscript "ie" stands for the "ideal evolution" , IIq is the initial field-line mapping, 
while Ft and F_t are tangential boundary fiows, forward and backward in time, respectively 
(see Figure [2]). This means that for any footpoint P, its conjugate footpoint P = n5.'^(P) can 
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generally be found by following the lower path of the diagram 

P ^ P 

i ... T (11) 

Po ^ Po 

In other words, instead of tracing the field line that connects P and P at time t (see Figure [2]), 
one can first trace the trajectory of P backward in time to find its initial prototype Pq, then 
tracing from it the field line of the initial configuration find its conjugate prototype Pq, and 
finally tracing the trajectory from Pq forward in time find P. Such a three-step calculation 
of the field-line mapping allows one to avoid the use of magnetic field data except for the 
initial data. Thus, whenever the tangential boundary flows and the initial magnetic field are 
known, the squashing factor may be calculated for any time moment without computing the 
ideal MHD evolution of configuration itself. 




Fig. 2. — Field-line mapping IT^ is expressed at any time t in terms of the initial field- 
line mapping IIo and the tangential forward and backward boundary flows Ft and F_t, 
respectively [see equation ( fTOl) ]. This mapping occurs whenever the frozen-in law is valid 
for the evolution of the configuration. The field lines and trajectories of the footpoints are 
shown by solid and dashed lines, respectively. 



Th is is very important if one takes into account the results of num erical MHD simu- 
lations (ILongcope fc Strauss! 1 1994 I Aulanier et al.ll2005l : iTitov et al.ll2008l ) that demonstrate 
a good correlation between regions with a highly distorted field-line mapping and strong 
current concentrations. In the light of this, we conjecture that equation ffTOj) must provide 
a powerful tool for predicting favorable sites of current-layer formation in magnetic config- 
urations whose evolution is subject to the frozen-in law. Its first successful application to 
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a particular configuration with two twisted magnetic flux spots (iTitov et al.ll2008l ) strongly 
supports this conjecture. 

In coordinate notations, the fleld-line mapping deflned by equation (fTUj) is described as 




{w\w'), (12) 



where the coordinate functions and their values are denoted by the same letters but different 
cases — the upper-case letters represent the functions, while the low-case letters represent 
their values; the functions and coordinates referring to the initial time moment are marked 
with the subscript 0. Thus, for example, the flow F^t niaps (m^, u^) to (mq = Uq{u^,u^), = 
Uq{u^,u'^)) and so on. The second step in this composition requires solving of the problem 
described by equations ([7])-([9]), with the proviso that this problem is set for the initial 
conflguration, whose boundary coordinates are (mq,Mq) and (wq,Wq). The flrst and third 
steps of composition f[T^ require tracing of the footpoint trajectories, which are determined 
from the tangential components of the boundary velocity fleld Vh, assumed to be known. 
Speciflcally, one needs to solve the following two initial value problems. For the flrst step, 
the system 

(fo,fo')|,., = (14) 
should be integrated backward in time from t to 0, while, for the third step, the system 

"^^^ 9''(l^-^^ . = 1,2;. = 1,2; (15) 



dt \dw 

(W\W%^^ = {wlwlr, ' (16) 

should be integrated forward in time from to t. In both these systems, summation over 
the repeating index s is assumed; the contravariant metric [g^^] = refers to points of the 
launch and target boundaries for equations ( fT3l) and ( |T5l) . respectively. The double-asterisk 
symbol (in the superscript) in equation (fT6|) indicates that the point {wq,Wq) is the image 
of {u^,u'^) that is obtained at the flrst two steps of composition (fT2|) . In other words, the 
double asterisks provide a double pull back of this point to the launch boundary, where Q is 
evaluated. This convention will be used hereafter for other values as well, with the extension 
to a triple pull back denoted by triple asterisks, respectively. 



It is clear that the flrst and third steps in the procedure described above are no more 
difficult than the second step. In principle, making consecutively these three steps at least 
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three times in the small neighborhood of a given launch point, one can evaluate numerically 
the Jacobian matrix ([1]) and the respective value of Q [see eqs. ([2])-([6])]. In practice, however, 
the divergence of the field lines at QSLs of the initial configuration can be so high that it may 
significantly affect the accuracy of such an evaluation. A better control of the accuracy in this 
procedure is obtained by using the fact that the Jacobian matrix Die of the composite "ideal" 
mapping (fTOj) is the product of the Jacobian matrices of the individual mappings, namely 



(17) 



Ae = 


'dW'' 




duQ 








M" 











Thus, the calculation of D is reduced to the calculations of these three matrices at the points 
that "lie on the path" of the composite mapping. As stated above, each of these calculations 
is similar to one another, and none of them requires knowledge of the magnetic field, except 
at the initial moment. 

It is clear from this consideration that the respective squashing factor can be ob- 
tained now by modifying equations ([2])-(jl]) as follows: 

tr {DlG*^D^^ G-i) 



|det {DIG*** D,,G-^)\^'^' 



The determinant of the product of the individual matrices in this expression can be rewritten 
as a product of the corresponding individual determinants. Using also conservation of the 
magnetic flux in the initial configuration, we finally obtain 

Bll tr(DTG***AeG'-i) 



^0 9 
9o9* 



nO 



B. 



det(M** M- 



We used here the fact that det 



'gll {B^Q^/g^), where 5*o and B^q are the normal 
components of the initial field at the corresponding conjugate footpoints, with = det Gq 
and Qq* = det Gq* at such points. 



2.2. Nonideal evolution of the magnetic field 

The breaking of the frozen-in law causes a change of the initial magnetic connection 
of the boundary plasma elements. In other words, such elements experience, in this case, a 
slippage relative to their ideal MHD mapping. This intuitive notion can be formalized by 
introducing two types of mapping, which relate two time moments and t in the evolution 
of a given magnetic configuration. The first of them, denoted by St and expressed as 



St = F^t oUtoFto Ho, 



(19) 
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maps the boundary onto itself at the initial time 0. We will call this the slip-forth mapping 
because the slippage of footpoints occurs forward in time. The second mapping describes a 
corresponding slippage backward in time and it is expressed as 

S_t = Ft o Ho o F_t o Ut (20) 

and called the slip-back mapping. This expression also maps the boundary onto itself but 
at the final time t, so that the slip-back mapping is an antipode to the slip-forth mapping. 
A pointwise definition of these mappings is depicted in Figure [3] for the case of a closed 
magnetic configuration. The picture can be easily modified, however, to represent more 
general configurations comprising both closed and open field lines by simply enclosing the 
coronal volume under study with the help of an additional upper boundary surface. 




Fig. 3. — Slip- forth (a) and slip-back (b) mappings [see equations (|T9l) and (|20i) . respectively] 
act on the footpoints of the field lines in such a way that each footpoint P and its image P 
are generally different due to a breakdown of the frozen-in condition. 
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Comparing equations (120|) and ffTOj) . we obtain 

= o n^, (21) 

which means that the shp-back mapping is just a composite of the "real" and "ideal" field-line 
mappings at time t. Similarly, equation f|T9l) can be rewritten as 

St = njfono, (22) 

where 

n|f = F^toUtoFt (23) 

is that initial field-line mapping whose ideal evolution driven by the fiow Ft would produce 
at time t the final field-line mapping lit. Thus, both slip mappings identify the differences in 
the field-line connectivities of the "real" and "ideal" configurations. The slip-forth mapping 
St, however, does this for the initial observer, while the slip-back mapping S-t does the same 
for the final observer. 

This consideration also shows that the footpoint slippage that will occur or has occurred 
within time t is fully determined by the magnetic fields at the initial and final moments and 
by the respective displacements of the plasma elements due to tangential boundary fiows. 
Such a slippage is caused by both magnetic reconnection and resistive plasma diffusion, 
whose contributions may generally be comparable in value. Y et it is clear that ma gnetic 



reconnection is more localized in space than resistive diffusion (ISchindler et al.lll988l ). One 
can expect therefore that these two processes are discriminated not by the magnitude of 
the slippage itself but rather by its spatial gradient. Note also that the introduced slip 
mappings are very similar to the field-line mapping in the sense that all of them are just 
two-dimensional mappings of the boundary on itself. Thinking by analogy, we can assume 
then that the squashing factor is a major characteristic not only for a field-line mapping but 
also for slip mappings. This assumption is fully confirmed by our further analysis, which 
shows that in the case of slip mappings the squashing factor makes it possible to identify 
the reconnecting magnetic fiux tubes in any evolving configuration, even if its evolution 
involves a substantial resistive diffusion. In particular, section H] presents an example of 
magnetic evolution in a current layer with a small nonideal region inside. It shows that 
the transition from magnetic diffusion to reconnection is not abrupt but continuous. The 
reconnecting flux tubes here pass through, and in the vicinity of, the nonideal region and 
they are distinguished, as expected, by high values of the squashing factor that are related 
to a locally enhanced current density. 
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To derive the required expressions, let us first represent the shp mappings in coordinate 
notations. For the shp- forth mapping defined by equation f|T9l) . we have 



St : iuQ,Uo) 



no Ft Ot F~ 

1 2\ (W-O'^o) / 1 2\ (w\w^) , 1 2\ / 1 2\ C^O'^o) / 1 2\ 



[u ,u 



(24) 



so that its Jacobian matrix is determined by the foHowing product of individual Jacobian 
matrices: 



D 



sf 



idun 






* 


'dUP' 






* 


'dWf 






duP _ 




dw'i 











(25) 



This composition implies that the respective squashing factor Qsf can be derived from suit- 
ably modified equations (l2])-(j4j) as follows: 



Q 



sf 



tr (D^GfAfGo- 
det (DjfGf AfGo' 



1/2- 



Rewriting the denominator in this expression as a product of individual determinants of 
the entering matrices and using conservation of a magnetic flux in the initial and final 
configurations, we finally obtain 



Q 



sf 



9*9o 



1/2 



9**9o 



BnBnO 



tr (DjfGf D^fGo^) 



(26) 



det {M-^***M*) 

Similarly, the slip-back mapping defined by equation ( l20l) is represented in coordinates 



by 



Ft 



S^t-{U,U) > [W ,W ) ^ {Wo,Wq) ^ {Uo,Uq) > (m ,M ), 



(27) 



SO that its Jacobian matrix is determined by the following product of individual Jacobian 
matrices: 









* 































(28) 
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This composition implies that the respective squashing factor Qsh can be obtained from 
suitably modified equations ([2])-(IlI) as 



-1 



Q 



sh 



tr (^DlG**DshG 



|det {DlG^** AbG-i) 



1/2- 



Rewriting the denominator in this expression as a product of individual determinants of 
the entering matrices and using conservation of a magnetic flux in the initial and final 
configurations, we finally obtain 



Q 



sb 



9o9' 
9o*9* 



1/2 



R * R* 



tr (D^G-AbG-^ 
det (M***M-i*) 



(29) 



By construction, the Qsf factor is defined on the boundary of the initial magnetic con- 
figuration by its subsequent evolution within time t, while the Qsh factor is defined on the 
boundary of the final magnetic configuration by its preceding evolution within time t. Both 
slip-squashing factors depend only on the two Jacobian matrices of the field-line mapping. 
Do and D, corresponding, respectively, to the initial and final moments. 

To better understand the meaning of these characteristics, consider a situation where 
the variation of the footpoint positions due to nonideal processes in the volume is much larger 
than that due to the surface flows at the boundary. From the physical point of view, this is 
actually the most interesting situation, because it corresponds to the case where magnetic 
reconnection going on in the volume is so quick on the time scale under study that the 
boundary motion of plasma can be neglected. The whole paths of the slip mappings are 
formed in this case only by field lines of the initial and final configurations, while Ft and F_t 
in (fT9!) - (!20l) are approximately identity mappings characterized by unit Jacobian matrices. 
Taking into account these observations, one can approximate equations ( !26l) and (!29l) by 

Qsf ~ 
Qsh ~ 

which show that, in general, if some of the matrix elements in Dq or D are large in absolute 
value, Qs{ or Qsh should be large as well. Such matrices, however, are related to the Q factor 
of the initial or final magnetic field in a similar way. Thus, Qst or Qsh should be large at the 
mapping paths [see eqs. fl2^ and fl271) ] that include the field hues belonging to the initial or 
final QSLs of the evolving configuration. 



9**9*0^ 


r 




9*9oy 




BnBnO 


9o*9*^ 


r 


TO** r?* 


9o9**y 




BnoBn 



tr [D'^D'^* G** D*Do Gq ^) , 
ti [D'^D'^* G** D*DG-^) , 



(30) 
(31) 
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By analogy with the QSLs, one can expect that the flux tubes assembled from such field 
lines have a layer- like structure. We will call these flux tubes reconnection fronts (RFs). Their 
determination is of great importance for analyzing magnetic reconnection in complicated 
3D configurations. Such configurations, however, require special consideration, and so the 
discussion of them is postponed to section 12.31 We note here only that both shp mappings 
have two RFs, one of which, called the instantaneous RF, corresponds to the present moment 
(0 or t, respectively, for the slip- forth or slip-back mappings). Equations (!30|) and (13T!) suggest 
that such a reconnection front should approximately coincide with the set of QSLs existing 
at a present moment (see an example in section [3]). However, the situation is different for 
the second RF, which is called the future RF or past RF, respectively, for the slip-forth or 
slip-back mapping. The future RF consists of those initial flux tubes whose footprints will 
turn after reconnection into the footprints of the final QSLs. The past RF consists of those 
final flux tubes whose footprints are formed by reconnection from the footprints of the initial 
QSLs. These definitions of RFs are slightly modified if the advection of footpoints is essential 
in the evolution under study. 



2.3. To-be-reconnected and reconnected magnetic fluxes 



It might seem surprising at first glance that our description of magnetic reconnection 
does not invoke an electric field at all, even though the componen t of the electric field 



parallel to a magiietic f ield is at the cornerstone of the GMR theory (jSchindler et al.lll988 



Hesse fc SchindleiUlOSsh . The physical ground for this omission of the electric field becomes 
clear if one recalls that its parallel component appears generally in very strong current 
layers (CLs), in particular, when the corresponding resistive voltage drop is not negligible. 
However, QSLs are likely to be favorable sites for the formation of strong CLs under very 



2003 



general conditions (Longcope fc StraussI 1994 : Priest fc Demoulin 1995; Titov et al 
Galsgaard et al.l l2003l : lAulanier et al.l l2005l ) . This makes it possible to put QSLs at the 



core of general reconnection theory by reformulating it in purely geometrical terms. Such 
a reformulation has the advantage of being able to directly relate the reconnection process 
with an accompanying variation of the structure of evolving magnetic configurations. More 
details illustrating these general statements are given in sections [3] and |H 

What is particularly important is that our approach enables us to identify reconnecting 
flux tubes in a given magnetic evolution as well as to estimate the reconnected fluxes by 
tracing the footprints of the RFs. To describe this more precisely, let us introduce first the 
lines going along the ridges of the Q distributions by denoting such lines as QSS and RF, 
respectively, for the field-line and slip mappings. QSS abbreviates here a quasi-separatrix 
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surface, which is a magnetic surface passing through the ridge of the respective Q distribution 
at the bo undary. QSSs pla y the role of the middle surfaces for QSLs [see some examples of 



QSSs in (ITitov et al.ll2008l )]. so that QSS lines serve as midlines for the footprints of QSLs. 



Similarly, the RF lines serve as midlines for the footprints of RFs. 

Consider first the simplest case, where the advection of the footpoints is negligible in 
comparison with their displacements due to reconnection, which means that both Ft and F_t 
are approximately identity mappings in composition fl24l) . Assume in addition that there are 
QSLs at both time moments and that their footprints are sufficiently separated in distance. 
This implies that the squashing of the St mapping on its composite path passing through 
a given footprint of the QSL is mainly provided by the field-line mapping corresponding to 
this QSL. By noticing also that Hq ^ = Hq, we obtain then from relation (12^ that 

RFo ~ QSSj, (32) 
RF± ^ no( QSSf ), (33) 

where the subscripts or t correspond, respectively, to instantaneous or future RFs, while 
the superscripts + or — denote the positive or negative magnetic polarity, to which the RF- 
or QSS line belongs. 

These equalities establish a direct relationship between RFs and QSLs and provide a 
simplified algorithm for determining RFs in the described situation. As discussed in section 
[H QSSs and boundary surfaces do not necessarily bound closed volumes in the corona, and so 
RFq- and RF ^ lines might not form closed contours at the boundary. The latter may happen, 
however, if QSSs, for example, turn out to be NP SSs. In this case, the closed contours, 
formed by RFq - and RF^ lines, outline in the initial configuration the footprints of the fiux 
tubes that will reconnect within time t. We will call them to-he-reconnected magnetic fiux 
tubes. In a more general situation with open RF lines, the footprints of to-be-reconnected 
fiux tubes are the boundary areas that are swept by RF^ lines in the initial configuration 
with time running from to t. The integral of the normal component of the initial magnetic 
field over such areas determines the magnetic fiuxes that will be reconnected within time t. 

Thus, the RF lines play a crucial role in determining the reconnecting magnetic fiuxes — 
the more accurately we determine the geometry and kinematics of the RF lines, the better 
we can quantify the reconnecting fiux tubes for a given magnetic evolution. Therefore, it is 
important to understand also what happens if the reconnection is not fast enough to neglect 
the advection of footpoints. To understand this question, let us turn again to equation (l24l) 
and its more compact form (l22l) -(l23l). from where one can see that the location of RF q lines 
should be given by equation (|32l) if the paths of St [see eq. (1241) ]. on which the field-line 
mappings IIq or IIq are mostly distorted, are different. The latter is probably valid if the 
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time interval t under study is not too small and so the footprints of QSLs of the Ho or IIq 
mappings do not overlap each other. 

Yet the situation with the RF^ lines is more complicated, even if the latter is true. The 
answer in this case depends on whether or not the boundary flows Ft and F_t significantly 
interlock or unlock the field lines in the volume of the final magnetic configuration. This, in 
turn, depends on whether or not these flows contain substantial shearing motions applied in 
transverse directions to the conjugate footprints of QSLs. Because a twisting pair of such 
shearing motions is able to seriously modify the Iljf mapping compared to the Ut mapping. 
In fact, this kind of advection can even destroy the original QSL or, in contrast, i t can create 



a new one, depending on the sign of the twist and the resulting values of shear (ITitov et al. 



20031 : iGalsgaard et al.ll2003l : iTitov et al.ll2008l ). If such an interlocking or unlocking effect of 



the boundary flows is small (but the footpoint advection is still not negligible), the RF^ lines 
can be approximated by 

RFi^ ^ nooF_i(QSSJ). (34) 



Otherwise, one needs to use a more accurate approximation, namely, 

RFf ^ no(QSS;^), (35) 

where QSSj denotes the QSS lines of the IIq mapping defined by equation (l23i) . Thus, the 
sequence of equations fl33|) . flM|) and fl35l) provides more and more accurate approximations of 
the RFj lines of the slip-forth mapping with respect to its distortion by footpoint advection. 

Using equations fl7r|) . fl?I]) . and fl2Ul) . a similar sequence of approximations can be found 
for the RFs of the slip-back mapping. In the case of a negligible footpoint advection, one 
obtains that the RF lines of the instantaneous and past RFs, respectively, are given by 

RF^^ ^ QSS f, (36) 
BEi ^ n,(QSSJ). (37) 

If footpoint advection is not negligible but capable of producing only a moderate distortion 
of the slip-back mapping, the RF hues of the past RFs should be determined by 

RFo^ ^ n,oF,( QSS^ ). (38) 

The presence of substantial interlocking or unlocking shears in the boundary motion requires, 
however, even a more accurate approximation, namely, 

RF± ^ n,(QSSj), (39) 
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where QSS denotes the QSS hnes of the Ilf mapping in the corresponding polarities. With 
time running backward, the RF ^ hnes of the instantaneous RF wiU remain practically un- 
changed [see equation (|36l)]. while the RFq lines of the past RF will change and move on the 
boundary by sweeping the footprints of the flux tubes that have been reconnected within 
the time t under study. 

The evolving topology of these RF lines and footprints may be rather non-trivial, since 
the reconnection process may generally occur at different moments and places in the con- 
figuration by causing changes in its structure. Such an intimate link of reconnection to the 
magnetic structure is expressed by the above approximate relationships between the RF and 
QSS lines. The accuracy of these relationships can always be proved by using exact expres- 
sions (1261) and (J29l) for the slip-squashing factors. But even using these exact expressions, 
the reconnected or to-be-reconnected magnetic fluxes may be estimated only approximately. 
The reason for this is that, in general, there is no strict difference between reconnection and 
magnetic diffusion, because the diffusion is always a part of the reconnection process. In 
our formulation, this fact manifests itself in an inherent uncertainty of the end points of 
the unclosed RF lines. Such end points can be identified by some threshold value Qthrcsh of 
the slip-squashing factors. Physically, this value is the degree of squashing of the elemen- 
tal flux tubes that start at the end points of the RF lines and pass through the border of 
the reconnection region. Since the location of this border is defined up to a fraction of the 
thickness of the reconnecting current layer, the threshold values Qthresh and the end points 
of the unclosed RF lines also cannot be defined exactly. 

From the physical point of view, the most accurate method for estimating them seems 
to rely on the voltage drop at the initial and final moments along the magnetic field lines that 
belong to the paths of the slip mappings starting at the respective RF lines. Analyzing the 
distribution of this voltage drop along the RF lines would give us the most precise criterion 
for the end points of these lines — they would be simply the points where the voltage drop 
decreases to a fixed fraction of its nearest maximum at the RF line. However, such an 
approach requires additional computations involving additional data on electric fields in a 
given magnetic evolution, which seems to be hardly justifiable from the practical point of 
view. The final result obtained by this electric field method will still not be unique but 
dependent on the threshold value of the voltage drop that is selected. In this situation, it 
seems to be more practical to choose some large value of Qthresh for identifying the end points 
of the RF lines by keeping in mind that the final result depends on the chosen Qthresh- It 
should be emphasized that such an ambiguity of our method appears only in the case of the 
unclosed parts of the RF lines — their closed parts provide an unambiguous description of 
the reconnection process (see an example in section [3]). It is interesting that the closed RF 
lines seem to appear less often than the open ones. 
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There may also be the cases, where the transverse displacements of the RF lines are 
negligible, so that the RF lines evolve mainly by changing their length (see an example 
in section H]). The reconnecting magnetic fluxes in such cases should be determined via 
the corresponding RF footprints defined at the boundary as areas with the slip-squashing 
factors larger than Qthresh- Thus, this generalized definition of reconnecting flux tubes also 
does depend on Qthresh, which reflects an inherent property of the reconnection process in 
general 3D configurations. 



3. EXAMPLE 1. FORMATION AND RISE OF A FLUX ROPE 



Let us illustrate now our theory of reconnection by considering a particular example of 
magnetic fiel d evolution. For th is purpose, we choose the evolving configuration analyzed 
previously by lHesse et al.l (120051 ) with the help of the GMR theory, whose major element is 
a nonvanishing voltage drop along the reconnecting magnetic field lines. The consideration 
of this example will significantly facilitate a comparison between the two theories. 

In a Cartesian system of coordinates {x,y,z), the corresponding components of the 
magnetic field are given by 





= -1 


By 


= 0.2 




= X, 



2' 



(40) 

(41) 
(42) 



where Ly and Lz are len gth-scale pa r amete rs, and e{t) is a monotonically growing function 
of time t. In contrast to lHesse et al.l (120051 ) . the variables y and z are swapped, so that the 
plane z = corresponds to the photosphere. Also the power degree 2 is recovered in the 
denominator of the second term in B^ to correct a misprint in the original formula. Only 
the Bx component evolves with time via the function e{t) > 0, whose time derivative e is 
assumed here to be strictly positive. This implies the presence of an electric field with a 
nonvanishing component 



e{t)z 



(43) 



(i + yyiDii + zyLi) 

such that Faraday's law for the evolving B^ is satisfied. 

By construction. By is nonvanishing and constant, so that the magnetic field (H0l) - (H2l) 
has no null points. It is easy to prove also that, for the assumed 6{t) > 0, B^ < at the 
photospheric polarity inversion line defined hj x = z = 0. Therefore, the transverse field 
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at this line is alw ays directed from the positive to negative polarity, which, in turn, means 



( iTitov et al.lll993l ) that this configuration has no coronal field lines touching the photosphere. 
Thus, under these assumptions, the topology of magnetic field lines remains trivial during 
the evolution of the configuration. 

Nevertheless, its magnetic structure changes dramatically, which can be immediately 
seen from the evolution of the transverse magnetic field B trans = {Bx,Bz) in the middle 
plane y = 0. Note that the field lines of -Btrans in this plane are contours of the function 

+ + (44) 

Figure m shows several such field lines at two different instances corresponding to e = 8 = 
and e = 13 = 62- They clearly demonstrate the formation of a magnetic island from the 
initial arcade-like configuration within the respective time interval. This 2D plot provides, of 
course, only a first glimpse of the evolution of the magnetic structure, which is in fact fully 
three-dimensional, since the vector field B depends on all three coordinates. This suggests, 
however, that the minimum points of might help to characterize the evolving magnetic 
structure. In this approach, the formation of the magnetic island in the 2D plane can be 
viewed as a bifurcation of the cusp-type minimum into elliptic and hyperbolic minima. This 
point of view is confirmed by a closer inspection of equations fHOl) - fHTl) . 



Indeed, one can see from them that the minima of B^ coincide with the nulls of Bx, 
whose locations are determined explicitly by 



± 



2e 



1 + yyLi \8{i + yym 



- 1 



1/2 



1/2 
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The upper and lower signs at the inner square root of this expression correspond to the loci 
of elliptic and hyperbolic minima, respectively, with heights zq and zx- The heights have 
real values only if £ > 8 and \y\ < l/max, where 

e \ 1/2 

o - 1 ) • (46) 



Ly 



As e increases from ei = 8, a single minimum point of cusp-type (yc = 0, zc = V^L^) 
bifurcates into two curves of elliptic and hyperbolic minima (see Figure [5]). It will be clear 
soon that some of these minima play a key role in the nonideal evolution of this configuration. 



In a highly conducting plasma, the perpendicular component of an electric field E_i 

.yyij - Byl 



Ey{y — ByB / B"^) mainly sustains the advection of plasma elements with a drift velocity 



V = E± X B/B . In our example, this component vanishes at the photosphere {z = 0) 




Fig. 4. — Field lines of the transverse magnetic field (5^., B^) in the middle plane y = at 
two different instances of the evolving magnetic configuration defined by equations P0|) - fj42|) . 



y/Ly 

Fig. 5. — Loci of minima of at four different times in the plane a; = 0. The sohd and 
dotted hnes, respectively, represent the minima of hyperbolic and elliptic types which emerge 
with growing e from the initial cusp point, present at £ = 8. 
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together with the total electric field Ey [see equation (H3|l ]. which implies that the pho- 
tospheric plasma elements remain at rest during the evolution of the configuration. It is 
of particular importance that, in the volume, there is a nonvanishing parallel electric field 
E\\ = EyByB/B'^, whose presence in the configuration provides a slippage of the footpoints 
with respe ct to the stat i onary photospheric plasma elements. The latter was clearly demon- 



strated by iHesse et al.l (120051 ). who have also shown that such a slippage is most intensive 
at the magnetic field lines with the highest values of the voltage drop, called by them pseu- 
dopotential. The latter is defined by 

S = - y E|| d/, (47) 

where / is an arc length of the magnetic field lines along which E\\ is integrated between their 
end points, so that S is a function of these poin ts. For convenienc e, we have reproduced the 



photospheric distribution of S, the same as in (IHesse et al.ll2005l ). except that, in our case, 
E\\ was integrated along the entire magnetic field lines without a cutoff by the lateral and 
top boundaries of numerical box. This was done to eliminate the artificial distortion of the S 
distribution due to such a boundary effect. The resulting S distributions at e = ei and e = 62 
are shown, respectively, in panels (a) and (b) of Figure El The second distribution differs 



only slightly from that calculated by iHesse et al.l (120051 ) . which means that the mentioned 



boundary effect is not large for the chosen size of the numerical box. 

For the same values of e, panels (c) and (d) in Figure El present the photospheric distri- 
butions of the Q factor in a logarithmic scale with the gray shading saturated at log Q = 4.0. 
Depending on the value of e, the actual maxima of logQ in these distributions are two-three 
times larger than the saturation level. So the largest values of Q are concentrated here 
in extremely narrow ridges tracing the footprints of the corresponding QSSs. Comparison 
between the Q and S distributions in Figure El [panels (c) and (d) vs. (a) and (b)] shows 
that the footprints of the QSSs pass exactly along the lines of the steepest gradient of S, or, 
in other words, along the lines where the spatial variation of the footpoint slippage is the 
highest. This indicates that the respective QSLs are indeed instantaneous RFs, as inferred 
in section [23] from the general argument. Further confirmation of this fact is obtained by 
considering the corresponding distributions of the slip-squashing factors. 

Before starting such a consideration, let us look first at the structure of the magnetic 
field lines in the vicinity of the QSSs. Figure [7^ shows a set of the field lines at e = ei with 
the launch points chosen to be very close to the cusp minimum point of B^. It is clearly 
seen there that the footpoints of these field lines are distributed exactly along the ridges 
of the Q distribution. So such field lines approximate the QSS by showing that it has a 
wedge-like shape with the sharpest corner at the cusp minimum point. Panels (b) and (c) 
in Figure [7] depict a set of field lines at e = 62 with the launch points near the loci of the 
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Fig. 6. — Photospheric distributions of pscudopotential S (normalized to E{t)) [(a), (b)], 
squashing factor Q [(c), (d)] at e = 8 and e = 13 (left and right panels, respectively) 
and slip-forth and slip-back squashing factors Qsf and Qsb, respectively, between these two 
instants. 
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hyperbolic minima of B'^ defined by equation (HSl) . The footpoints of these field lines are also 
distributed exactly along the ridges of the corresponding Q distribution. So such field lines 
approximate the QSS by demonstrating that it consists of two parts attached to each other. 
The lower part has a wedge-like shape with a sharp corner along the loci of the hyperbolic 
minima. The upper part is a tube-like surface enclosing a twisted flux rope that lies on the 
top of this wedge. It starts at the photosphere as a narrow, very fiat, and almost untwisted 
flux tube spreading along the slopes of the wedge and blows up near the wedge corner into 
a very large, elongated, and highly twisted plasmoid. The cross section of the resulting QSS 
by the middle plane y = coincides with the separatrix line of the transverse magnetic field 
shown in the right panel of Figure |H Thus, the QSLs in our example are due to the presence 
of the hyperbolic or cusp minima of — this fact will be explained in section [5] from the 
general point of view. 

In the strict sense, the oval-like parts of the QSS lines (Fig. [6li) are not closed but 
open lines — the ridges of the Q distribution here approach very close to one another when 
making their turns but they never completely merge. We believe that t his wa s also the 



case in the twisted configuration analyzed previously by iDemoulin et al.l (1l996bl ) with the 
help of the norm N rather than Q. Such a behavior of QSLs in both cases has a simple 
explanation: due to a relatively large twist of the rope, there are field lines with a different 
number of windings (larger than one) in the neighborhood of the hyperbolic magnetic minima 
in the corona. A relative proximity of such field lines to these minima is translated into 
the apparent closeness of the respective QSS lines at the photosphere due to an enormous 
squashing that the elemental flux tubes experience on their way from the minima to the 
photosphere. From the observational point of view, however, this seems to be a atypical 
situation, since the observed flaring bright loops rarely show multiple windings. Therefore, 
the QSL footprints in realistic magnetic c onfigu r ations should look like open rather than 



closed curves, as actually demonstrated by iTitovl (l2007al ) for a flux rope configuration with 
a realistic amount of twist. 

Consider now the distributions of the slip-forth and slip-back squashing factors calcu- 
lated numerically for the variation of e in the interval [ei,£2]- These calculations are signif- 
icantly simplified due to the above-mentioned line-tying condition E\^^q = that implies 
the absence of plasma flows at the photosphere. As discussed in the previous section, the 
calculation of Qst or Qsb in such a case requires knowledge only of B at the initial and final 
moments, which can be obtained here from equations (H0l) - (l42l) . The resulting distributions 
of Qsi and Qsh are shown in panels (e) and (f), respectively, in Figure O The inner or outer 
pairs of the ridges of the Qgf distribution trace the instantaneous or future RF lines, respec- 
tively; and vice versa: the outer or inner pairs of the ridges of the Qsh distribution trace the 
instantaneous or past RF lines, respectively. Comparing panels (e) and (f ) vs. (c) and (d) in 
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Fig. 7. — Magnetic field lines passing in the vicinity of the minimum points of the cusp 
[(a), e = 8] and hyperbolic types [(b) and (c), e = 13] and the respective photospheric Q 
distributions. To-be-reconnected [(d)-(f)] and reconnected [(g)-(i)] magnetic flux tubes at 
the initial {e = 8) and final {e = 13) moments, respectively; their launch footprints are shaded 
in orange and superimposed on the corresponding photosperic Qsi and Qsh distributions. The 
respective orientation of the coordinate axes are shown by red (x), green (y) and blue (z) 
arrows. 
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Figure [6] shows that the instantaneous RF hnes practically coincide with the instantaneous 
QSS lines — exactly as previously inferred from the general consideration [see equations fl32l) 
and fl5B]) ]. In contrast, the future and past RF lines are noticeably distorted in comparison 
with the corresponding QSS lines, which is also in accord with the inferred general equations 
(!33|) and (1371) . Such a distortion is sufficiently smooth, that it causes the oval- like parts 
of the corresponding RF footprints to appear nearly as closed as their corresponding QSL 
footprints. 

It is useful also to consider the behavior of the RF lines in response to the variation of 
£2 or El for the slip-forth or slip-back mappings, respectively. In both cases, the footprints of 
the instantaneous RF lines remain practically unchanged under such a variation by recurring 
the shape of the corresponding QSS lines shown in panels (c) and (d) in Figure El On the 
contrary, the future and past RF lines change their locations, shapes, and even topology. 

Let £1 = 8 be kept constant and 62 be increasing in value from 8 to 13. Then the future 
RF line in each of the magnetic polarity will split and move outward from the respective 
instantaneous RF line by gradually increasing its length and bending itself in the middle part. 
In addition to these transformations, two small and very flat ovals emerge in the future RF 
lines by stretching themselves out from the conjugate footpoints of the field line that passes 
in the volume through the cusp minimum point. Such ovals enclose the footprints of the 
initial fiux tubes that will reconnect at the moment 62 and form the above-mentioned flux 
rope — these flux tubes are shown in Figure [TJi. The areas swept out by the major parts of 
the future RF lines determine the footprints of the flux tubes that will reconnect and form 
the "flare loops" below the flux rope. One can see from Figures [7^ and[7f that, as expected, 
their conjugate footprints also include the magnetic flux that will belong to the flux rope. 
Finally, the areas that are swept by the ovals in the exterior determine the footprints of the 
flux tubes that will pass through the flux rope formed during this evolution. 

Let 62 = 13 be kept constant and ei be decreasing in value from 13 to 8. Then the 
past RF line in each of the magnetic polarity will split and move inward from the respective 
instantaneous RF line by gradually decreasing its length and contracting the ovals that 
enclose the footprints of the flux rope (see Figure [7^). At the moment ei = 8, the ovals 
will completely shrink to the conjugate footpoints of the fleld line that passes in the volume 
through the cusp minimum point. The areas swept out by the major parts of the past RF 
lines determine the footprints of the reconnected flux tubes manifesting themselves as "flare 
loops" formed below the flux rope [panels (h) and (i) in Figure [7j . This kind of flux tube 
may have been directly observed in X-rays as coronal loops slipping along their respective 



QSL footprints (lAulanier et al.l 120071 ). Finally, the areas that are swept by the ovals in the 



exterior determine the footprints of the flux tubes that have passed through the flux rope 
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formed during this evolution. 

Thus, this example clearly demonstrates that the slip-squashing factors allow one to 
identify reconnecting magnetic fluxes in a given magnetic evolution with an unprecedented 
level of detail. 



4. EXAMPLE 2. RECONNECTION IN A CURRENT-LAYER PATCH 

Let us consider now an alternate magnetic evolution, where, in contrast to Example 1, 
the footprints of the initial and final QSLs overlap each other. This causes the resulting RFs 
to behave in a less predictable way than before, and it also provides some new insights into 
the significance of the slip-squashing factors. 

As the initial configuration, let us take the force-free magnetic field depending only on 
the X coordinate as 

Bq = sm6{x)y + cos6{x)z, (48) 
e{x) = eotanh{x/lsh)- (49) 

The corresponding magnetic field lines are straight and tilted toward the z-axis at the angle 
6{x), which changes smoothly with x from —6q to 6q to form a sheared structure. Most of 
the magnetic shear, and its associated electric current density, 

are concentrated in a slab-like layer of half thickness Igh- Let this configuration be bounded 
from the top and bottom by ideally conducting plates located at 2; = ±L = ±1. Then 
using the straightness of the magnetic field lines and equations ([I])-([5]), it is not difficult to 
calculate the respective squashing factor 



cosh'*(x//sh) cos^6'(x) 



Q = 2+ I / (51) 



Figure El^a) shows the profiles of Q at different s-nd other fixed parameters by demon- 
strating that at sufficiently small Ish and large 9q our current layer becomes a QSL, in which 
both Q and current density increase together with the total shear angle 29q. Since |-Bo| = 1 
everywhere in the configuration, the appearance of such a QSL is not related to the presence 
of any magnetic minima in contrast to Example 1. Instead, this particular QSL is the result 
of a large magnetic shear present inside a thin layer with a strong current density. 
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Fig. 8. — Profiles of Q across the initial slab configuration (e = 0, /gh = 0.05, = 0.025 
and Ly = Lz = 0.3) at 6o = kn/lS, with integers k = 1, ... ,6 (a). The distributions of Q 
at z = —1 (b) and z = (c), the distributions of — S/e (d) and Qsh (e) at 2; = — 1 in the 
perturbed slab configuration with e = 0.006 and 6q = tt/S (the remaining parameters are the 
same); the contours of B'^ — 0.91 + 0.01 A;, A; = 0, . . . , 9, are superimposed on the respective 
Q distribution in panel (c). Panel (f) shows for the same parameters the gray-shaded areas 
in the plane z — —1 where logQsb > 2 for four different times corresponding to £ = 0.001 A; , 
A; = 2, . . . , 6, with a stepwise decrease of shading from the initial to the final times. 
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If the latter reaches a certain threshold value defined through other plasma parameters, a 
nonideal plasma effect may come into play. For example, the onset of some plasma instability 
might be triggered. This would, in turn, cause the appearance of an electric field with a 
nonvanishing component along Bq inside the layer. One can assume also that in reality some 
of the plasma parameters may have nonuniform distributions and so the threshold conditions 
would be reached first at some spot in the layer. This would, in turn, lead to a localization of 
the parallel electric field in the spot. We can model such a process in our case by introducing 
into the configuration the following electric field 

E = e{t)f{x,y,z)z. (52) 

The time and space variation of E is described here by the time derivative of some function 
e{t) and the function 

fix, y, z) = [x'/Ll + y^Ll + cosh(z/L,)] , (53) 

respectively. Taking the curl of E and integrating it over time, we obtain in accord with 
Faraday's law the resulting magnetic field 

B = Bo + Ezx V/. (54) 

This is, of course, a very approximate model of the respective real process, since E is 
prescribed here in an ad hoc way rather than derived from the appropriate plasma dynamics. 
However, the model is accurate enough to qualitatively demonstrate the difference between 
magnetic reconnection and diffusion in a simple, but generic, case. 

The field-line equation does not admit an analytical solution for the derived magnetic 
field B, so the corresponding squashing factors have been calculated numerically for the 
following set of parameters: Oq = tt/S, Ish = 0.05, = 0.025, Ly = = 0.3, and e = 0.001 k 
with integers A; = 1 ... 6. In all these cases, the assumed electric field is localized in a patch- 
like region that has the shape of a very fiat ellipsoid compressed in the x direction. Figure 
m^b) shows an example of the Q distribution in the plane z = —1 at e = 0.006, when the 
transverse gradients of the perturbed and unperturbed y components of the magnetic field 
become comparable in value at the center of the patch. The respective Q distribution at the 
top plate 2; = +1 is obtained from the one below by its mirror refiection about the y-axis. 
These distributions show that two well-pronounced QSLs adjacent to one another eventually 
emerge inside the initial QSL. Such an evolution of the QSLs is very different from the one 
studied in Example 1. 

By definition, Q is invariant to the direction of mapping along the field lines, so its 
footpoint value can be prescribed by any other point of the respective field line, thereby 
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extending the boundary Q distribution into the volume. This allows one to study the 3D 
structure of the QSLs via the distribution of Q at different cross sections. Figure[8](c) displays 
one such distribution in the plane ^ = at £ = 0.006 together with the respective set of 
contours of B"^ = const in this plane. One can see from this figure that the positions of 
local maxima of Q and minima of in the plane z = are relatively close to one another, 
although not identical. A noticeable shift between them is probably the result of the relative 
shallowness of the minima of B^, whose values are just 10% lower in comparison with the 
background, while the characteristic sizes of the minima are of the order of Igh- In other 
words, the appearance of the new QSLs seems to be related again to the formation of the 
minima of B^ in the configuration, as was the case in Example 1. The only difference is that 
such a relation in Example 1 is stronger because of a much larger depth of the corresponding 
minima. 

For the same set of parameters. Figure [8](d) shows the distribution of the pseudo- 
potential S [defined by equation fHTl) ] in the plane z = —1, which suggests that a substantial 
change of magnetic connectivity must occur in our configuration along the diagonal of the 
computational rectangular area. This is confirmed by Figure IH](e) that presents the respec- 
tive distribution of the slip-back squashing factor Q^h for the time period, where e increases 
from to 0.006. Since the assumed electric field drops exponentially with growing \z\ toward 
the boundary plates [see equations fl52|) - fl53l) ]. the Qsb distribution is calculated here under 
the assumption that the line-tying condition is strictly fulfilled at the plates. In other words, 
a possible small footpoint advection is neglected in calculating the respective slip-back map- 
ping, which is composed of the final and initial field-line mappings only, as described in 
section 12. 2[ The Qsb so obtained demonstrates that a major change of connectivity occurs 
along a narrow strip that is aligned with both the respective new QSL footprints and the 
area with large values of The latter is in good agree ment with the ana lysis of a similar 



nonideal evolution in the so-called X-line configuration (IHesse et al.ll2005l ). The indicated 
narrow strip can be identified here only with the footprint of the instantaneous RF, since the 
past RF is absent at the initial moment (simply because the new QSLs just start to emerge 
at that moment). 

To see how this RF footprint evolves, we determine it as the area at the boundary 
plate, where Q^h > Qthresh = 10^, for several pairs of e, whose initial values equal 0. This 
value is used same for all the pairs, while the final values vary within the interval [0, 0.006]. 
These RF footprints, represented in Figure [8]^f) by gray shadings of different intensity, form 
a clearly nested structure such that the area that is swept out by the final time approaches 
asymptotically to the area of the instantaneous RF footprint at this time. In terms of 
the corresponding RF lines, this means that they grow in length by remaining nearly self- 
aligned and not by moving in the transverse direction (in contrast to Example 1). Therefore, 
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in order to determine the reconnected flux tubes in situations like this, one should extend 
our previous definition of the reconnected boundary area by replacing the sweeping RF 
lines in this definition with the sweeping RF footprints. Since both of them depend on the 
chosen value Qthresh, such an extension does not increase the ambiguity in determining the 
reconnected magnetic fiuxes. It does, however, make this definition more general. 




Fig. 9. — Magnetic field lines in the slab configuration aX e = 0.006 with other parameters, 
the same as in Figure [SI The launch points are located in the planes z = —1 (a) and z = +1 
(b) along the contours logQsb = 2.0. The respective orientation of the coordinate axes are 
shown by red (x), green (y), and blue (z) arrows. The length scale in the x-direction is 
stretched ten times. 



The reconnected fiux tubes as defined above are shown in Figure |9] with the help of the 
magnetic field lines that are launched from the contours logQsb = 2.0 at the boundary plates 
z = ±1. These fiux tubes intersect but do not coincide with one another, which is evident 
from the locations of the target footpoints being scattered partly outside of the RF footprint 
at the other boundary plate. Nevertheless, such footpoints also outline a very elongated 
oval that is slightly turned with respect to the corresponding RF footprint around their 
common center (a; = 0, y = 0). This is a manifestation of the characteristic plasma rotation 
in magnetic fiux tubes passing through a noni deal region — aii effect that was previously 
described in the framework of the GMR theory (jSchindler 2006 ). 
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Our second example suggests that the characteristic feature that distinguishes recon- 
nection from diffusion is RFs with narrow ribbon-hke footprints. Such footprints are charac- 
terized by two length scales, one of which is much smaller than the other. Displacement of 
plasma elements in the nonideal region across the corresponding RFs causes a large change 
in magnetic connectivity along the RF footprints. In other words, the high values of slip- 
squashing factors at the flux tubes that pass through and in the vicinity of a nonideal region 
is the major signature of reconnection compared to diffusion. The transition from one pro- 
cess to the other is not abrupt but continuous both in time and in parameter space, and only 
in this sense it is just a matter of semantics of how one refers to this process — reconnection 
or diffusion. Generally, the appearance of such nonideal regions in configurations is closely 
related to locally enhanced current density, which corresponds to large magnetic shear, large 
values of the squashing factor, and hence large values of the slip- squashing factors. Because 
of this causal chain, it seems preferable to refer to the process in the nonideal regions as 
reconnection rather than magnetic diffusion. 

Note also that strong outflow jets of plasma from the nonideal regions have been a 
feature of traditional 2D reconnection models. We may conjecture that perhaps this is also a 
property of many 3D reconnection reg imes, although perhap s some types of 3D reconnection 



do not possess rapid outflow jets [e.g. ISchindler et al.l (Il988| )]. The accurate proof of such a 



conjecture must rely on a self-consistent analysis of the full system of MHD equations, which 
is beyond the scope of the present work. However, some indications of its validity can be seen 
from our present consideration as well. In particular, equations flS5]) and flMl) show that the 
considered nonideal process causes an X-type perturbation of the transverse magnetic field. 
Therefore, the Lorentz force due to interaction of this perturbation with the initial strong 
current layer will tend to accelerate plasma out of the nonideal region and to produce outflow 
jets. Thus, the thinner the current layer, the stronger the acceleration, which corresponds, in 
turn, to the narrower RFs developing in the configuration, as demonstrated in this section. 



5. MAGNETIC MINIMUM POINTS 

Determining QSLs and RFs via computation of the corresponding squashing factors is 
conceptually simple but technically difficult. However, the particular examples of magnetic 
evolution considered in sections [3] and H] strongly suggest that there is a close relationship 
between QSLs, RFs, and magnetic minimum points. If the presence of such a minimum 
point is a sufficient condition for the existence of the respective QSL, this might significantly 
simplify determining QSLs in a given configuration. Finding magnetic minima is techni- 
cally as simple as finding magnetic null points, so locating such minima would significantly 
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facilitate the localization of QSLs and RFs. This motivates us to study the local properties 
of magnetic minima in the most general terms. 



Bq ^ 0, the system of coordinates can always be 



Note first that, for any point with 
chosen in such a way that B^j. = By = 0, while B^ = Bq. If the point is a magnetic minimum, 
then Vi?^ = V-B^ = or simply V-B^ = in such a system of coordinates. Taking also 
into account that V ■ -B = 0, we obtain the following general expression for the matrix of 
gradients of magnetic field at its minimum point: 

/ d^Bx dyBx dzBx \ 
dxBy —dxBx dzBy 


Rotating the system of coordinates about the z-axis does not change the components of B 
at the minimum point, but it reduces B at a certain rotation angle to the form 
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where all matrix elements are generally different from those in equation ([55] 
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Thus, in 

this rotated system of coordinates, we obtain the canonical linearized form for a magnetic- 
minimum field 



By 

B. 



ax — by + cz, 
bx — ay + dz, 
Bo. 



(57) 
(58) 
(59) 



Using these expressions, one can easily find from equation ([7]) the corresponding magnetic 
field lines. Depending on the values of parameters a and b, there are only three different types 
of minimum: hyperbolic, cusp, and elliptic minima, exactly as in the example considered 
in section [31 Assuming that the field lines start in the plane z = 0, so that {x,y, z)\^^q = 
{xq, yo, 0), we have for all these types of minimum 
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which simply represents linear dependence of z on the differential flux tube volume r. For 
hyperbolic minima, which exist when \a\ > \b\, the field lines are given by 



X 



, , , sinh UT cosh ur 

Xq + {axo - byo) ■ h cBq- 



sinh/^r cosh /ir 

yo + [bxo - ayo) h dBo — 

/i /x^ 



1 / 7 7\ n sinh ur — ur 

- + {ac-bd)Bo ^ , ^ , 61 

1 sinh ur — ur , , 

- + (6c - ad)Bo ^„ ^ , (62) 
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1 /9 

where fx = (a^ — IP') . In the hmit /i — yielding \a\ = one can retrieve from here the 
field-hne structure at cusp minima described by 



X = xo + a{xo - yo)r + cBo— + a{c- d)Bo—, (63) 

2 o 

y = yo + a{xo-yo)T + dBo— + a{c-d)Bo — . (64) 

/ 6 

The same expressions can also be retrieved in the limit uj ^ from the elliptic minima, 
which exist when \a\ < \b\ with neighboring field lines determined by 

X = xo + (axe - by^f^ + cB,'-^^^ + (ac - 6rf)i^o """^"' ^\ (65) 

,, , sin cur ,^ cos cur— 1 , , , ,^ sin cur — cur 
y = yo + {bxo-ayo) \- dBo \-{ad-bc)Bo , (66) 

CU CU"' (jJ-^ 

where u = {P — a^)^''^ . 

The above analytical expressions show that the highest divergence of field lines in re- 
sponse to the variation of Xq or yQ is achieved in the vicinity of hyperbolic magnetic minima 
[see equations (IM]) and (1^ ]. The corresponding gradients of x and y with respect to Xq or 
yo grow exponentially with z, and the rate of this growth increases with decreasing Bq. In 
the case of cusp minima, however, such a growth is only linear in z, while for elliptic minima 
it is even saturated via the sine-function. This provides a strong general argument in favor 
of the idea that the hyperbolic and cusp minima of a magnetic field are responsible for the 
appearance of QSLs and corresponding RFs in evolving magnetic configurations. 

Although both of the examples considered in sections [3] and H] possess some symmetry, 
we do not think that this is a crucial factor for the emergence of QSLs and RFs. The really 
important parameters that control the thickness of the emerging QSLs should be the depth 
of the respective hyperbolic magnetic minima and their characteristic sizes compared to the 
global length scale of the configuration. In this respect, the magnetic null points are just a 
degenerate case of magnetic minima, whose associated QSLs are transformed into genuine 
separatrix surfaces. More detailed studies are certainly needed in the future to quantify the 
relation between QSLs and magnetic minima. 



6. SUMMARY 

We have developed a general theory for the evolution of magnetic connections between 
boundary plasma elements in an arbitrary 3D magnetic configuration with an arbitrarily 
shaped boundary. The theory assumes that at each instant the magnetic field is known 
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everywhere throughout the volume domain including the boundary, where the tangential 
flow of plasma elements is also known. 

The major quantity that is use d for characterizing the magn etic connectivity and its 
evolution is the squashing factor Q (ITitov et al.lll999l : lTitovll2007al ). After being applied to 
the magnetic field-line mapping, the Q factor allows one to identify in a given configuration 
all separatrix and quasi-separatrix surfaces, as well as their hybrids. We have demonstrated 
first how the Q factor is determined for any ideal evolution of the magnetic configuration 
in terms of the initial magnetic field and tangential boundary flows without reference to 
the magnetic field at other moments in time. This result is itself of great importance for 
the further development of the theory of current-layer formations in magnetic configurations 
with a highly conducting plasma via its motion at the boundary. 

In order to describe the nonideal evolution of magnetic connectivity, we have introduced 
the so-called "slip mapping" , which characterizes the change of magnetic connections between 
the boundary plasma elements due to violation of the frozen-in law in the plasma volume. 
This mapping can be defined either forward or backward in time to give correspondingly slip- 
forth or slip-back mappings, which operate on the initial or final configurations, respectively, 
within a chosen time interval. By calculating the squashing factor for these new mappings, 
we have defined new quantities called slip-forth and slip-back squashing factors, Qsf and 
Qsh, respectively. The large values of these quantities identify the reconnecting magnetic 
flux tubes, called reconnection fronts (RFs), at the initial and final time instances. The 
boundary areas that are swept by RFs within a chosen time interval correspond to the 
footprints of reconnecting magnetic flux tubes. In this way, the Qsi factor determines in 
the initial configuration the to-be-reconnected magnetic flux tubes, while the Qsb factor 
determines in the final configuration the reconnected magnetic flux tubes. The dependence 
of the properties of the RFs on tangential boundary flows is also analyzed. 



The general results are illustrated with two specific examples b ased on evo. 
netic configurations 



The first of these was considered earlier by iHesse et al .■ _ 
the framework of the general magnetic reconnection (GMR) theory (jSchindler et al 



vmgmag- 
J2005h in 



1988 



Hesse &: Schindlerlll988l : iHesse et al.lll99ll : iHesse fc Birnlll993l ). Revisiting this example, we 
have demonstrated that our approach is consistent with GMR theory but complements it in 
several important aspects. The second example describes a magnetic evolution that is driven 
by a thin current layer with a field- line voltage drop at a small patch inside it. This example 
enables us to show that the localization of such a nonideal process in a small region may 
indeed cause the structural changes that are characteristic for reconnection rather than for 
magnetic diffusion. The latter example, in particular, demonstrates the formation of mag- 
netic flux tubes with narrow, ribbon-like footprints, which is a typical feature of reconnection 
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process. 

In comparison with GMR theory, our approach handles regions with small or large 
amounts of reconnected flux on an equal footing. This is because the slip- squashing factors, 
in contrast to the voltage drop of GMR, are dimensionless geometrical quantities measur- 
ing only a relative spatial rate of magnetic slippage at the boundary. For the same reason, 
the slip-squashing factors discriminate between regions of reconnection and simple diffusion 
much better than the voltage drop method. Moreover, since the slip mappings are composed 
from the field-line mappings and both types of mappings are characterized by the same quan- 
tity, the analysis of reconnection based on our theory is intimately related to the analysis 
of magnetic structure. As we have demonstrated here, the latter provides intriguing oppor- 
tunities for analyzing magnetic reconnection in 3D configurations with an unprecedented 
level of detail. The use of slip-squashing factors opens new perspectives for the comparing 
reconnecting flux tubes predicted by the theoretical models with the different morphological 
features observed in the solar flares and coronal mass ejections [ e.g. brightening coronal 



loops, chromospheric kernels and ribbons, see the recent review by iDemouliru (120071 )]. 



Our analysis of these examples also reveals the special role of hyperbolic and cusp 
minimum points in 3D magnetic field configurations. The elemental magnetic flux tubes 
passing in the vicinity of such points may become extremely squashed, which implies the 
formation of the respective QSLs. We have extended these examples by studying the field- 
line structure in the neighborhood of generic magnetic minima and confirmed this particular 
result. Hence, similar to magnetic null points, hyperbolic and cusp minima provide an 
important localization of plausible reconnection sites. More detailed investigations in this 
direction promise to lead to even more interesting results. 

It should be emphasized that, similar to GMR theory, our approach is applicable to the 
analysis of reconnection in magnetic configurations described not only by MHD, but also 
by any other plasma model. The algorithm of calculating slip-squashing factors admits a 
straightforward implementation and it is highly parallelizable. Our experience shows that, 
even without parallelization, the squashing factor may be computed for known magnetic 
field data with a modest demand on computing resources. Thus, the whole approach is 
practical enough to pursue the implementation of the proposed algorithms for calculating 
slip-squashing factors for the solar corona. The solar physics group at SAIC is presently 
working on this project. 
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